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Abstract 

The  Kohn-Hulthen  variational  prjrriple  is  expressed  in  terms 
of  a  scattering  operator  formalism.  It  is  shown  that  in  terms  of 
this  formalism  the  principle  can  be  expressed  more  abstractly  and 
concisely  than  heretofore. 
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1.  Introduction 

A  ▼arlatlonsl  principle  has  "been  devised  "by  Hulthen  for  phase  shifts  in 
scattering  problems  (ref.  l).   The  principle  has  heen  modified  and  its  appli- 
es tions  hare  been  extended  by  Kohn  (ref.  2),   Both  of  these  authors  have  set 
up  the  variational  principle  explicitly  in  terms  of  scattering  phases.  How- 
ever, it  is  simpler  to  set  up  the  variational  principle  in  terms  of  the  scatter- 
ing operator  formalism.  Furthermore  in  terms  of  this  formalism  the  principle 
is  more  abstract  and  more  general  than  in  terms  of  the  phase-shift  representa- 
tion.  In  order  to  discuss  this  more  abstract  formallsmj  we  shall  sumnsrize 
the  properties  of  the  scattering  operator,  using  essentially  the  concepts  of 
Miller  (ref.  3),  Triedrichs  (ref.  ^),  and  Lip^^man  and  Schwinger  (ref.  5). 

2,  The  Scattering  Operator 

In  scattering  problems  the  Hamiltonian  H  of  the  system  is  broken  up  into 
two  perts,  one  of  which  Is  called  the  unperturbed  ftimiltonian  H  and  the  other 

of  which  is  celled  the  perturbation  eV.  Thus 

(1)  H  «  H  +  eV. 

o 

H  is  80  chosen  as  to  have  a  continuous  spectrum. 

I'et  us  denote  the  eigenvalues  of  H.  by  E.   Generally  the  spectrum  of  H 

0  0 

is  degenerate.   Wp  shall  therefore  introduce  accessory  dynamical  variables  which 
commute  with  H  and  which  together  with  H  form  a  complete  set  of  commuting 
variables.  We  shall  denote  the  eigenvalues  of  these  accessory  variables  col- 
lectively by  a.  We  shall  then  write  the  eigenfunctions  of  H  corresponding  to 

0 

the  eigenvalue  E  and  a  as  «  (E,a),  where 

0 

(2)  H  (B  (E,a)  =  E«)  (E,a.) 

(?)  («  (E.a),  0.  (E«,a«))  =  6(E-B  •)6(a,a«); 

0       o 
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the  left-hand  side  of  ('^)  means  the  Hermltian  Inner  product,  and  6(a,a*)  is 
a  generalized  8-function  defined  "by 

(i!)      y  r(a)6(a.,c')da  =  r(a')  if  a»  Is  in  »2 

=  0  otherwise, 

where  J   is  the  ran/re  of  integration.   Integration  over  those  variables  in- 
cluded in  a  which  lie  in  the  discrete  spectran  are  to  "be  replaced  by  siuinnationa, 
It  might  "be  useful  to  indicate  the  nature  of  the  ▼arlahles  a  in  order  that 
they  njBiy  he  thought  of  in  more  concrete  form.   If  one  considers  the  scattering 
of  a  particle  by  a  perturbing  potential,  one  may  take  H  to  be  the  kinetic 
energy  of  the  particle  and  cT  may  be  considered  as  the  scattering  potential. 
In  this  case  a  might  be  chosen  to  be  the  polar  angles  which  specify  the  direc- 
tion of  the  momentum.   One  might  also  consider  the  scattering  of  an  electron 
by  an  atom.  H  would  then  be  the  sum  of  the  kinetic  energy  operator  of  the 
electron  and  the  Haailtonlan  of  the  atom.   In  this  case  a  would  consist  of  the 
polar  angles  which  give  the  direction  of  the  electron's  momentum,  together  with 
the  vsriables  which  denote  the  state  of  the  atom.  For  en  example  of  a  problem 
of  this  kind  see  Moses  (ref.  6), 

We  shall  denote  an  arbitrary  state  J  in  the  H  representation  by  f  (l!,a) 
and  call  it  the  "representer"  in  that  representation;  we  write 

(5)  I<T>    ^^^'"^  "  («o(lB.a,).  i  ) 
or 

(5a)  '§_mJjf{-E,a.)   (i)^(E,a)dadB  . 

An  operator  K  in  the  H  representation  will  be  written  as  K^;  it  is  defined  by 

(6)  K  5<-~>K^  f(E,a). 
In  particular 

(7)  ^  =  B 

as  follows  immediately  from  (2)  and  (5a),   Sometimes  it  is  convenient  to  write 
&  a?  an  integral  operator  with  the  kernel  or  "matrix"  i:^(E,a|E',a' ): 
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(8)  K^f(i:,rT,)=^K^(E,a|E',a')f(E',a')da'dJ!«. 
where 

(9)  E^CE.all'.a')  »  (•^(35.<t.),Ku.^(E  •  .a,t)). 

▼or  purpoBefl  of  the  scptterin^  operator  we  are  interested  In  ei^enfimctlons 

of  ^  which  correspond  to  eigen-ralnes  lying  in  the  same  range  and  having  the  same 

degeneracy  as  the  eigenvalues  of  H  ,  We  shall  denote  the  eigenvalues  of  H  hy 

o 

7,     To  descrlhe  snch  eigenfonctlons  of  H  we  shell  have  to  introduce  operator* 

which  commute  with  H  and  which,  together  with  H,  form  a  complete  set  of  commuting 
variables  collectively  hy  P.  The  values  of  3  have  the  same  range  as  a.   It  is  not 
generally  possible  to  give  P  a  direct  physical  meaning,  hut  it  Is  not  necessary 
to  do  so  for  our  purposes. 

The  eigenfunctions  of  H  are  not  unique.  We  shall  he  interested  in  two 
different  sets  of  eigenfimctions  which  we  shall  denote  hy  at  (5',P)  and  o(~(P,3); 
these  we  shall  now  characterize. 

Let  us  denote  the  •"  in  the  H  representation  hy  u  : 

(10)  n^(E.alP.P)  =  (u)^(E,a).  «r(F.p)), 

The  eigenfijnctlons  w    are  characterized  hy  the  following  equations  for  u^  : 

(11)  n^(i:,a|l',P)=  8(E-r)6(a.P)+  €V_^(F-Il)T^(ll,a|ir,p), 
where  V^(t)   is  essentially  6^(0: 

(12)  V^(0  =     lim         -^    =  ±itT6(0  +  P  7  . 

The  letter  P  in  (l2)  means  thPt  in  integrations  over  integrands  involving  t 
as  a  factor,  the  Cauchy  principal  part  is  to  he  used.  Also 

(13)  T_^(l!,CT,|F.P)=yy  7^(l!,alE',a')u^(B',a'|P.p)da»dll'«  7^u^(l,a|F,P), 
where  the  operator  v  operates  on  u.  insofar  as  u.  are  functions  of  E,a, 
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The  Bcattering  operator  and  its  Inverse  can  "be   o'btalned  In  terms  of  the 

fanctions  u.   as  follows:  Consider  solutions  of  £(t)  of  the  Schroedinger 

equation 

il^)  i  ^  E(t)  =  H5(t)      (t  =  1  ). 

If  the  limits 

iH  t 

(15)  lim   e  °  I(t)  =  5+ 
t->too 

exist,  then  the  scattering  operator  S  is  defined  "by 

(16)  £+  =  S  I,  . 

If  we  consider  the  H  -representer  of  J^t): 

0 

l(t)  <^>    f  (E,a;t). 
we  see  that  two  possihle  solutions  to  the  Schroedinger  equation  can  he  written: 
d'')  f^(1'-.a;t)  =y^u^(lS.n|r,e)e"*'V(r,P)dPdT, 

where  g(F,P)  is  a  quadraticfllly  integrehle  function  which  was  chosen  to  corres- 
pond to  a  choice  of  initial  conditions  on  f^  as  seen  helow.     Friedrichs    (ref,  ^4) 
shows  that  f     is  characterized  hy 

(18)  lim         e^*f    (E,a;t)  =  e(E,a); 
t->  -  00 

(19)  lim         e^*f    (E,a;t)  »//  S^(E.a|E',a«)g(E'a')da«dIS'. 
t->  4-  00  ~  ^"^ 

where 

(20)  s''(E.a|E«,a')  =  6(E-E')6(a,a')-2TTic6(E-E')T_(E'.alE«,a'). 

We  can  identify  equations   (18)  and   (19)  «8  "being  equation  (15)  in  the 
H  -reuresentatlon.     In  pnrtlculer,   if  we  identify  g(E,a)  as   the  H^-representer 
of  5      and  //s    (E,a|E«,a»)g(E ',a' )da'dE'  as  the  H^-representer  of  5^.,*®  B©e 
that   s''(l!,rTlE«,a»)   is  Just  the  kernel  of  the  Integral  operator  which  gives  the 
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ocptterln^  operator  S  In  the  H  -representation. 

0 

In  a  slmllfir  way  it  can  "be   shown  that 

(21)      lim    e^*f^(15.a;t)  -  g(E.a) 
t->  +  oo 

(?2)      Urn    e^V(l^.n.;t)  =  //  S^(E,alE',a«)g(S',a«)da'dJl« 
t->-oo     *        '/'' 

where 

(23)    S^(I!,a|i:«,a«)=  ^CE-E  •  )6  (a.a')+  2nic8(E-E  •)T^(E«  .aJB ',o»). 

It  is  clear  that  S  (E,a|E',a')  is  Just  the  kernel  of  the  integral  operator  which 
gives  the  inverse  scattering  operator  S~  in  the  H  -representation. 

In  order  to  obtain  the  scattering  operator  we  need  only  find  T  (E,a|E,a*)» 
We  shall  fonnalate  the  variational  principle  for  this  function.   It  should  "be 
noted  that  T  (E,alE,a')  is  the  amplitude  of  the  scattered  wave  in  the  x- 
representation.   It  is  easier  to  work  in  the  H  -representation  "because  the 
Green's  functions  involved  are  simply  V  -f-unctionsj  in  the  x- representation 
they  ere  considera"bly  more  complicated. 


''.  The  Vwristional  Expression 

The  Kohn-Hulthen  vnriptionnl  principle  is  "based  upon  variations  of  the 
following  functional  A(^  ): 

(2ll)  A(^^)  =  I((^^)  *  €^(E,P'l5'.e). 

where 

(25       I((7.)  =  lim   //v.(E,a|E',P')(H^-E)v  (E,a|F,P)dadE 
*    E  '->  E 

(?6)   v^(E.o,lF,p)=6(E-E)6(a.,8)+€V^(r-E)^_^(E,a|E,$) 

=  u^(E,a|E,p)+€V_^(E-E)[^_^(E,a!E,p)-T_^(E.a|r,pa  . 


The  "bar  means  the  complex  conjugate, 

fLvAftA    as    tT»4Al    •^»^Tlct.^  nns    ^m«  '_ 


The  functions  J   ^   are  to  "be  regarded  aa  trial  functions  for  T  .  W© 


shall  prove  later  that  variations  of  A( /^)  which  result  from  variations  of 

± 


"X  .    a"bout  T     aVe  zero.      Hence  ACJ    )   is  an  Insensitive  functional  of  T 
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in  the  vicinity  of  s'  .  =  T  and  we  can  write 

(27)  M?^)     -  A(T^). 

where  J        approximates  T  .   It  will  also  "be  shown  that 

(28)  I(T^)  =  0 

and  that  hence  relation  (27)  can  he  written 

eT^(?,P«|?.P)  =  1(?^)*€   -IT  (F.P'If.P)  =  A(7_j_) 

or^  eqtii\ralently, 

(29)  T_(F.P«|F,P)  =iA(7^). 

Relation  (29)  is  the  form  of  the  Kohn-Holthen  principle  in  terms  of  the 
scattering  operator  fomelism.  We  use  trial  f-unctions  J  in  the  right-hand 
side  of  (29)  to  ohtaln  p  more  accurate  expression  for  T^(T,P '|F,P). 

It  should  he  noted  that  in  evaltiating  the  expression  (25)  for  l(^_j^) 
one  must  take  care  not  to  interchange  the  limit  and  the  integration  processes 
prematurely.   If  one  takes  the  limit  hefore  carrying  out  at  least  pert  of  the 
integration,  the  integrand  is  not  defined  hecause  one  has  the  product  of  two 
symholic  fimctions  ▼.  v  which  have  the  same  arguments  in  terms  of  the  varia- 
hles  S  and  I*.  Products  of  symholic  functions  having  the  same  argument  ere 
generally  not  defined;  however  if  the  symholic  functions  have  different  argu- 
ments  their  products  can  he  defined.  For  example,  (8(x))  is  not  defined,  hut 
6(x)6(y)  is. 

After  some  manipulation  it  is  possihle  to  hreak  up  the  integrand  into  a 
sum  in  which  terms  involving  products  of  symholic  functions  are  isolsted.  The 
integral  over  these  products  cpn  he  evaluated  and  the  limit  T'  — >  F  is  then 
taken,  for  the  remaining  terms  it  is  proper  to  interchange  the  integration 
«d  limiting  processes.  The  procedure  will  he  illustrated  in  the  proof  of  the 
variational  procedure  pnd  in  the  discussion  of  some  results. 
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i'.     Proof  of  the  TTohn-Hiilthen  Varle.tlonal  Procedtire 

We  stell  flrBt  prove  equBtlon  (28),     We  have  from  the  definition   (25) 
of  I(^^)  and  definition    (26)  for  v 

(30)  T(T   )  =  lim       ffnATi,r\T\e*)(^^)u   (E,a!F.e)dadI. 

However,  we  remember  that  in  the  H  -representation  u  (E,a|F,P)  la  an 
eigenfunction  of  H  and  n     is  the  operator  H,   Hence  Ho)~(F,P)=  ^^"(F.p)  becomes 

(31)  H^u_(E,alF,p)  =  Fu_(E.alF.P) 
or 

(32)  (^-F)u_(E,alr,P)  =  0. 

From  this  It  follows  that 
(28)  I(T^)  =  0. 

We  shall  now  show  that  variations  of  ^ij  .)  restiltin^  from  independent 

variations  of  J        a"hout  T     vanish. 

We  denote   this  variation  of  HT    )  hj  6A(cr    )L     «    ;  then 

±     ± 

(33)  6A(v^^)      =  A(T^+6T^)  .  A(T^). 

Halng  (?M)t    (25),  (26)  we  have  to  the  first  order  in  6T^ 

6A(s/^  )  =  lim  cE/7u.(E.alF«.3')(H^-F)V  (F-E)eT  (E,alF.p)dcadI5 
-   F»->F   •'•'  "^  ~       ~ 

(3i^)  +^V^(F«-E)6T_^(F,a|F«.P')(H^-F)u_(E.a|F,P)dadJl] 

4-  €  6T_(F,P'|P,P). 

The  second  term  in  the  right-hand  side  of  aqtiation  (3^)  vanishes  becatise  of 
equation    (32). 

I«t  us  designate  the  first  Integral  on  the  right-hand  side  of  (3**)  "by 

^'(7   ),  which  we  now  evaluate. 
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(35)  A'(T  )=  11m  e  //  5.  (E,a|y»,p' )(H^-5')v  (r-E)8T  (E.a|r,P)dadil. 

_  Y  I— >  f    JJ      "*"  "       ~ 

Since  H  is  a  Hermitian  operator  and  since  u.  Is  an  elgenfunction  of  H  in  the 
H  -representation, 

(36)  H^u^(E,alF'.e')  =  J'u^(E.a|r',p')} 
therefore  we  hare 

A«(T  )  =  lim   €  //  Ij:h^-F)u  (E,a|r',P'2Y  (PJ:)6T  (E.a|F,P)dadIl 
i   j>  t— >.  Y     -^J  ~  " 

^^"^^        =  lim   €  /7u^(E,aIr'.e«)(r'-E)V  (r-E)6T  (E,aIr.P)dadB. 

We  note  that 

(Et_T')v  (p_E)  =  QF'-E)-(F-E2*^  (J'-E) 

(T8)  =  (F«-E)V  (?_E)-(r-E)V  (F-E) 

=  (F«_E)V  (F.-E)-!. 

When  (38)  Is  used  in  (3*^)  w«  have 

A«(T  )=  lim  eC/7't[.(E,alF«,P«)(F'-E)V  (F-E)6T  (E,o,|r,p)dadB 

-  F«~>F   '''' 

^^^  -^u^(E.alF'.e')6T_(E,a|F,p)dadE;]. 

Now  from  (ll)  it  follows  that 

(^0)  (F-E)u^(I.a|r,p)  =  €T^(E,a|F.p). 

Hence 

A'(T  )=  lim  t\c   /7t  (E,a,|F«,P')V  (F-E)6T  (E,a|F,P)dadJi 

-  F  •— >  F 

(^1)  -y/'u^(F,a|F',p>)6T_(E,a|F,3)dadSl]. 

Note  that  in  the  first  of  the  integrals  we  no  longer  have  products  of  symbolic 
functions,  since  T  and  6T  are  well-behaved  functions.  Only  V  is  a  symbolic 
function.   In  the  second  integral  n  is  the  only  symbolic  function.  Hence  we 
can  now  Interchange  the  limit  and  integration  processes  in  both  integrals. 
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It  was  not.  possible  to  do  this  earlier.  Then  recognizing  that  V  (t)  =  V.(t), 
we  ohtain  from  (iil) 

(^2)  -^n^(E.alF.B')6T_(E.a|F.e)dad2l]. 

Wow  from  (ll)  it  Is  seen  that 

e\('-l!)5^(2.n|F.P')  =  n^(E,alF.p«  )-6 (E-r)6(a,P«). 
^ence 

A'(T^)  =  en^u^(^.^l''.P')6T_(E.a|r,e)dad2 

(^3)  -  6T_(F,P«|r.p)  -^u^(B,aIl!',P')6T_(E,a|F,p)dadEj 

=  -€  6T_(F,p»Ir.p). 

Finally,  since 

tk{J\)      =  A'(T^)  +  €6T_(F,P«|F,P). 

we  have  from  (^3) 

which  completes  the  proof. 

5.  Fxplicit  form  of  A(T ^) 

Although  the  form  of  ^(T^)   given  hy  equations  (2^)  through  (26)  is  useful 
for  proving  the  stPtlonsry  properties  of  this  functional,  it  is  not  very  con- 
venient for  explicit  calculation.   One  can  substitute  expressions  (25)  and  (26) 
in  (2^)  and  ohtain  the  following: 

f^(T^)   «  €/(F,et|ir,e) 

+  c^i^/(F,p'|l,a)v_(F-E)i^_(E,c».|F,p)dadJ: 

+  JJ   7^(F.alF,p«)V_(F-E)r'(F.a|F,P)dadE 
-  _^  :?^ ^  (E .a|F ,  p  •  )V_ (F-E  );r^  (E ,al F .  p)dad}y> 

+  €^y/y|/'r+(2,a|F,p')V_(r-E)7^(E.a|E»,a')7_(r^«). 

T   (E«,a«lr,p)dadEda'dE>, 
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^e  derivation  of  (^5)  is  straightforward,  though  one  mast,  as  usual,  take  care 
in  the  use  of  the  limit  proceBs  involved  in  the  definition  of  l(J    ), 

We  shall  now  Indicate  a  few  simple  results.  First  of  all  we  note  that 
T  (E,a|r,P)  can  he  calculated  by  a  perturbation  procedure.  Eqxjationa  (ll) 
throtigh  (13)  can  he  combined  to  form  an  integral  equation  for  either 
u  (E,o.|F,e)  or  T  (E.nly.P);  this  eqtiation  can  be  solved  by  interation.  First 
one  uses  the  approximation  T  5  0  and  inserts  this  in  equation  (ll)  for 
u  (E,«.|F,P).  With  this  value  of  u  ,  one  uses  (l?)  to  obtain  better  approxima- 
tions for  T  ;  this  procedure  can  be  repeated  to  obtain  any  desired  degree  of 
accuracy.  This  method  is  equivalent  to  solving  for  u  or  for  T  by  using  a 
power  series  erpanslon  in  e.  Fach  iteration  used  to  find  T_^  is  called  a 
Bom  approximation  for  T  ,  For  example  the  zeroeth  Bom  approximation  for 
T_^(E,a!F,$)  is  zero.  The  first,  v  (E,a!F,P),  is  obtained  by  using  the  first 

mm 

Bom  approxlmstlon  for  u  ,  namely  u_^(F,a|F,p)=  6(E-F)6(a,P),  in  (13). 

let  us  use  the  zeroeth  Born  approximation  for  T^  as  a  trial  function*'  ♦ 
That  is,  let  us  use  as  trial  functions  T^   =  0.  Then  we  have  from  (^5)  and  (29) 

(46)        T^(F,3»|5'.P)  =  v^(r,p«|F.p). 

This  is  Just  the  first  Bom  approximation  to  T_(F.P'|^.P)  that  would  be  obtained 
rising  thp  Iterstion  scheme  discussed  ^bove. 

Let  us  take  as  trial  functions  V   the  first  Bom  approximation  to  T  , 
namely 

'^+(E,a|F,P)  =  /(E.alF,P). 

Substituting  into  (^^5)  and  using  the  fact  thet  V  Is  a  Hermitian  operator,  we 
have 
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T_{?,fi«|r,e)  =  /(T.e'!T,P>+<^V^(y,3'i^.a)^.(y-S)  V^(B,a|F,e)dadJ 

This  Is  jvLBt   the  third  Born  approximation  to  T  (F,P'|?,e)  using  the  iteration 
procednre.  Kohn*  points  out  thet  if  one  uses  for  J      the  m-th  Bom  approxima- 
tion to  T_^  and  for  J      the  n-th  Bom  approximation  to  T  ,  then  the  variational 

expression  giree   T  (T.P'If.P)  accurately  to  the  (nr»-tt+l)  approximation.   It 

•hould  he  noted,  however,  that  the  variational  expression  in  no  way  can  he 

used  as  an  iterative  procedure  to  obtain  T  ,  since  the  Hulthen  procedure 

gives  an  approximate  function  to  T  (E,a|l',P)  only  for  E  =  F,  whereas  one  needs 

trial  functions  which  are  defined  for  general  values  of  E  and  I.for  an  iterative 
method. 


*  W.  Kohn,  Private  communication. 
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